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Abstract 



Extensive numerical simulation are reported for the structure and dynam- 
ics of large clusters on metal(lOO) surfaces. Different types of perimeter hop- 
ping processes makes center-of-mass of the cluster to follow a a random walk 
trajectory. Then, a diffusion coefficient D can be defined as lim D(t), with 

t—tOO 

D(t) = (d 2 ) /(At) and d the displacement of the center-of-mass. In the simula- 
tions, the dependence of the diffusion coefficient on those perimeter hopping 
processes can be analyzed in detail, since the relations between different rates 
for the processes are explicitly considered as parameters. 
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The different diffusion processes that take place on surfaces have, clearly, an important 
role in many technological areas. Diffusion of individual atoms and clusters has been studied 
for a long time with different experimental techniques and, more recently, using scanning 
tunneling microscopy (STM) |14]||. From a theoretical point of view, several lattice-gas 
kinetic Monte Carlo simulations |§ have addressed the question of the dependence of the 
cluster diffusion coefficients on the cluster size in square HJ^] or triangular lattice. Re- 
cently, the diffusion of metal clusters on metal surfaces has received extensive attention due 
to the fact that some experimental and theoretical studies have led to the expectation that 
only small, two-dimensional (2D) clusters were able to diffuse. The larger two-dimensional 
clusters, also observed on the surfaces, were not expected to diffuse. However, recent exper- 
imental evidence from STM studies @ have became available showing that very large 2D 
Ag clusters clusters (containing N = 10 2 — 10 3 atoms) are indeed able to diffuse on Ag(100) 
substrates. 

From the previous theoretical and experimental work, it is now clear that the movement 
of the cluster is a consequence of several atomic-scale processes, taking place at the periphery 
of the cluster, that makes the center-of-mass follow a random walk trajectory. For instance, 
for small clusters (N < 20) it has been proposed that the mechanism of diffusion is the 
short range motion of a single atom away from the periphery followed by a regrouping of the 
cluster around the peripheric vacancy || . Evidence of concerted gliding is also available || . 
On the order hand, the movement of clusters composed of a large amount of atoms has 
been considered a consequence of two main mechanisms taking place at the boundary of 
the cluster: periphery diffusion (PD) and 2D evaporation-condensation (EC). In the PD 
process several types of atomic motions along the periphery of the cluster are responsible 
for the displacement of the center-of-mass. However, the atoms do not leave the cluster 
while executing those movements. In the EC process the cluster is considered to be in 
quasi-equilibrium with a dilute 2D gas of atoms, diffusing very quickly on the metal surface 
surrounding the cluster. Both types of process are non-exclusive and can take place at the 
same time. Of course, it can be anticipated that the energy barriers for the PD mechanism 
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are much lower than the energy needed for an atom to leave the cluster and, in principle, 
the PD process can be expected to be primarily responsible for the movement of the cluster. 

At this point, the question of the dependence of the diffusion coefficient on the number 
of atoms in the cluster (N) becomes relevant. The diffusion coefficient D can be defined as 

\imD(t) = (d 2 )/(4t), (1) 

with d the displacement of the center-of-mass. 

It has been suggested that the value of the diffusion coefficient D behaves as D ~ 
N~ a . Different values for the exponent a have been proposed depending on which diffusion 
mechanism is considered to facilitate the movement of the cluster flT3"|,|Il |. In this sense, 



Monte Carlo simulations for cluster diffusion based on the PD mechanism are available, 



showing a = 1.5 to 2.0 [10,11]. However, these values of a lead to a strong variation of 



the diffusion coefficient D as a function of the number N, which is not consistent with the 



experimental STM data available 0]. On the other hand, in reference |1| it is claimed that 
for circular clusters the diffusion coefficient behaves as D ~ d~ n , where d is the diameter of 
the cluster and the integer n identifies different diffusion mechanisms. When the center of 
mass motion occurs by adatom diffusion along the periphery of the cluster n = 3, while n = 2 
or 1 when the cluster diffusion occurs by correlated or uncorrelated adatom evaporation and 
condensation, respectively. 

Then it is clear that a deeper understanding of the way in which 2D large clusters move 
can be obtained using numerical simulations. Following this idea, we set up Monte Carlo 
simulations for the kinetics in which each of the above described processes takes place at 
its own rate. Diffusion of atoms along the cluster perimeter proceeds at a rate r e and 
evaporation-condensation at a rate r c . Both are related by the detailed-balance relation 
i r c/ r e) — exp(— AE/kfyT), where AE is the energy difference between the two processes. 
The relation between the rates is kept as a parameter in the simulations. 

Our Monte Carlo simulation proceeds as follows. We start with a 'square' cluster of N 
atoms at the center of a 1500 x 1500 lattice substrate ||12||. At a randomly selected site from 



the periphery of the cluster two main events can take place. 

(a) With probability r c /r e , an evaporation or condensation event is selected. If the 
evaporation event is going to take place, an atom is detached from the cluster and goes into 
the surrounding two dimensional gas of fast, diffusing atoms. If the condensation event is 
going to take place, an atom is attached to a randomly selected place on the periphery of 
the cluster since it is considered to come from the surrounding two dimensional gas. In this 
way, the average number of atoms in the cluster in maintained equal to N. 

(b) With probability 1 — (r c /r e ), an atom is relaxed to the nearest-neighbor or next- 
nearest-neighbor empty site where it finds that more bonds are saturated. This is the 
mechanism that makes the cluster stay connected (i.e. not dissolve) in the r e 3> r c limit. 
We restrict ourselves within this limit letting r c /r e be less than 0.2. 

After a certain thermalization time, we start to record the position of the center-of-mass 
of the cluster with respect to the origin of coordinates considered to be at the center of the 
lattice representing the substrate. Time is measured in Monte Carlo steps per atom. 

The simulational results are shown in Figs 1 to 3. Results presented have been averaged 
over 100 independent runs and for 10 different cluster sizes containing from 121 to 961 
particles. 

In Fig. 1 is a representative plot of the overall behavior of the diffusion coefficient. It can 
be seen that, in agreement with the available experimental results JF], the simulations show a 
fast decay in the value of the diffusion coefficient at early times, for (r c /r e ) = 0.1 and 121 < 
iV < 961. After this fast decay the cluster reaches a diffusive random- walk-like behavior, 
demonstrated by the plateau at t — > oo. This figure confirms that the main characteristics 
of the movements of the cluster are well reproduced in the simulations. In particular, the 
very fast decay at early times has been associated with a certain back correlation effect, and 
it has been observed in previous cluster diffusion simulations where it was considered to be 



related to the cluster connectivity |1C] 



In Fig. 2 the diffusion coefficient D is plotted versus 1/N. From the plot it can be seen 
that D depends on iV as D = Do(r c /r e )N~ a with a = 1. The most notable characteristic 



is the dependence of the pre-factor Dq on the evaporation-relaxation relation r c /r e . The 
dependence can only be observed by a numerical simulation that allows full control of the 
relation. It can be seen that the cluster has a greater D for greater values of the relation 
r c /r e . This result is a little bit surprising. It tell us that, although the predominant processes 
at the periphery is the PD process, the EC process has a greater effect on the mobility of 
the cluster. Then, the simulations favor the EC mechanism as the primary one responsible 
for the movement of the cluster, in accordance with the available experimental data. 

From Fig. 3, it can be seen that the exponent a does not depend on the ratio r c /r e . By 



linear fitting the points of the plot, an average ITBJ value for the exponent a is obtained. This 



value is a = 1.092 ± 0.123 and differs from the value obtained from the experimental data, 
just reported in reference |7|. However, the experimental measurements of D as a function 
on N have been adjusted within a certain experimental error and cannot be considered to 
be conclusive. On the other hand the value of a ~ 1.0 is in accordance with the theoretical 
predicted value of reference |15j (n = 1) when the evaporation-condensation process is 
uncorrelated. 

Finally, in Fig 4 the behavior of the diffusion coefficient vs time (at early evolution 
times) is plotted in log- log plot for r e — 1.0 and three different values of r c . From this plot 
a very new characteristic of the movement is identified. The diffusion coefficient behaves as 
D(t) ~ t _/3 with (3 = 0.812 ± 0.021, showing a nontrivial scaling behavior at early times. In 
principle, it could be anticipated that the scaling behavior should be in close relation with 
the above mentioned back correlation effect. But, the explanation of this characteristic from 
a physical point of view needs more experimental and theoretical work. 

In conclusion. Results for the simulation of the diffusion of large Ag clusters over Ag(100) 
surfaces are presented. The simulations reproduce very well the main characteristic of the 
movement of the cluster according to the available experimental results. The dependence 
of the diffusion coefficient on the number of atoms in the cluster has been explored, as 
well as other features of the mechanisms by which the large clusters diffuses. Also, from 
the simulations new characteristics of the movement have been identified and suggest new 



experimental and theoretical work. A deeper understanding of the basic cluster diffusion 
mechanisms will have important consequences on the techniques that make possible the 
epitaxial growth. 
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FIGURES 

FIG. 1. The behavior of the diffusion coefficient along the total simulation time for N = 121 
(top curve) to N = 961 (bottom curve). The very fast decay at early time is in correspondence 
with the experimental results. 

FIG. 2. The diffusion coefficient D versus 1/N. Different slopes for different values of r c indicate 
that the prefactor Dq depends on the r c /r e relation. 

FIG. 3. Log-log plot of D versus 1/N. The average slope of the three plots gives the value for 
the exponent a = 1.092 ± 0.123. 

FIG. 4. Log-log plot of the diffusion coefficient vs time (at early evolution times) for r e = 1.0 
and different values of r c . The straight decay indicates an scaling behavior D(t) ~ t~~@, with 
P = 0.812 ± 0.021. 
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